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Abstract 

Classic and recent results for gauge effects on the properties of the normal- 
to-superconducting phase transition in bulk and thin film superconductors are re- 
viewed. Similar problems in the description of other natural systems (liquid crystals, 
quantum field theory, early universe) are also discussed. The relatively strong gauge 
effects on the fluctuations of the ordering field at low spatial dimensionality D and, 
in particular, in thin (quasi-2D) films are considered in details. A special attention 
is paid to the fluctuations of the gauge field. It is shown that the mechanism in 
which these gauge fluctuations affect on the order of the phase transition and other 
phase transition properties varies with the variation of the spatial dimensionality 
D. The problem for the experimental confirmation of the theoretical predictions 
about the order of the phase transitions in gauge systems is discussed. 



1 INTRODUCTION 

1.1 Superconductivity and gauge effects 

A remarkable example of gauge theory in condensed matter physics is the Ginzburg- 
Landau (GL) functional of superconductivity [T, "2]. The latter is invariant towards both 
global U{1) rotations of the order parameter field ip{x) and local gauge transformations 
(rotations) of the same field and the vector potential A[x) of the magnetic induction B. 
These two properties of gauge invariance define the global and local U{1) symmetries of the 
GL theory. In both cases the gauge group is a one-dimensional Abelian continuous group, 
f/(l). The spontaneous breaking of these global and local symmetries below the phase 
transition point allows for the appearance of the superconducting phases: the uniform 
Meissner phase and, under certain circumstances, the mixed (Abrikosov vortex [3]) phase. 
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While the Meissner phase ^ 0, (B) = 0] is a mere product of the breaking of the 
global f/(l) symmetry, the vortex phase [{ip{x)) ^ 0, {B{x)) ^ 0] , where the equilibrium 
field configurations of both the order parameter field iplx) and magnetic induction B{x) 
are spatially nonuniform, is a result of the spontaneous breaking of local U{1) symmetry. 
The local gauge is important also for the description of the magnetic field penetration 
in both Meissner and vortex phases. This penetration is described by an additional 
characteristic length, the London penetration length [H [2], which is different from zero 
only when the symmetry is broken, i.e., in the ordered phases, where the equilibrium value 
of is different from zero. The local U{1) symmetry is present when the superconductor 
is in external magnetic field, or, when this field is equal to zero but magnetic fluctuations 
exist and give sense of the vector potential A{x) as purely fluctuating field. 

The fluctuations of the relevant fields, ip{x) and A{x) in usual superconductors are small 
and can be neglected. This means that the GL functional can be investigated in the 
lowest order mean field (MF) approximation (alias, "tree approximation" [HE]). This is 
the usual way of treatment of the GL free energy, in particular, the investigation of the 
GL equations [HISIE]. Within the tree approximation, the phase transition from normal 
to superconducting state in zero magnetic field is of second order [U [2l Hj. For a long 
time this phase transition has been considered as one of the best examples of second order 
phase transitions, which has an excellent description within MF. But in 1974 Halperin, 
Lubensky and Ma (HLM) [6] showed that the magnetic fluctuations change the order of 
the superconducting phase transition in a zero external magnetic field {Hq = \Hq\ = 0)), 
i.e., the order of the phase transition from normal-to-uniform (Meissner) superconducting 
state at T^o = Tc{Hq = 0) (see, also Ref. [7]). Since then this fluctuation change of 
the order of normal-to-superconducting phase transition {HLM effect) has been under 
debate. After 1974 an overwhelming amount of theoretical research on this topic has been 
performed but up to now there is no complete consensus about the order of this phase 
transition. In this review we shall consider some aspects of this problem but we shall 
not be able to discuss or mention all relevant contributions. However, it is important to 
emphasize that the investigation of the properties of the superconducting phase transition 
is important for other areas of physics, too. Here we shall briefly enumerate and discuss 
several examples. 

1.2 Quantum field theory and other problems 

In elementary particle physics the gauge invariant theory similar to the GL theory of 
superconductivity is called Abelian Higgs Model [H [9l [TOl [TTl [121 [13] • The same global 
[/(I) and local U{1) gauge symmetries are present but the phenomenon of spontaneous 
breaking of symmetry occurs only for imaginary mass of the Higgs (ip—) field, which is the 
analog of the field describing the Cooper pair in a superconductor. In the absence of 
spontaneous symmetry breaking this model would describe an ordinary electrodynamics 
of charged scalars, but the situation becomes more interesting when the mentioned mass is 
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imaginary and the breaking of symmetry occurs. Now the symmetry breaking phenomena 
receive other names and another physical interpretation. The breaking of both global 
and local gauge symmetries ensures a mechanism of transformation of the two initial 
scalar fields - analogs of the components ip' and ip" of the complex field = %p' + iip", 
and two massless photon fields - analogs of the two independent components Aj of the 
vector potential A = {Aj;j = 1,2,3;V.A = 0} in a superconductor, to four massive 
particle fields: the so-called Higgs boson, which is analog of the spontaneous order I?/;! > 
in a superconductor, and three massive vector field components, i.e. a massive three 
dimensional vector field. The mass of this new vector field is proportional to the electric 
charge \q\ and magnitude IV'I of the Higgs field and, as a matter of fact, this is exactly 
the way, in which the London penetration length in a superconductor depends on the 
electron charge |e| and the modulus |'0| of the superconducting order parameter ip. Thus 
the spontaneous breaking of the local gauge symmetry leads to the formation of massive 
particles without spoiling the gauge invariance of the theory. This is called "the Higgs 
mechanism" . The latter plays a fundamental role in the unified theory of electromagnetic, 
weak and strong interactions (see, e.g., Ref. [T2| [13]). 

It is easy to see that there is something quite common between the phenomena of sponta- 
neous breaking of the continuous symmetries in superconductivity theory and in quantum 
field theory. The superconducting phase I'?/'! > is the exact analog of the Higgs boson in 
the Abelian-Higgs model, whereas the appearance of a massive vector field has its analog 
in the finite London penetration length, mentioned above. Of course, there is no obsta- 
cles in interpreting the phenomena of spontaneous symmetry breaking in quantum field 
theory as phase transitions by taking the Higgs mode mass as tuning parameter. The 
phase transition will occur at zero mass of the Higgs boson. 

A similar phenomenon of spontaneous breaking of both global U{1) and local gauge sym- 
metries is possible also within the scalar electrodynamics due to mass insertions from 
the radiation corrections, as shown in Ref. [H]. The radiation corrections, analog of the 
magnetic fluctuations in a superconductor, generate an imaginary mass to the initially 
massless scalar field in this theory and the latter becomes very like the Abelian-Higgs 
model. Here the symmetry breaking leads to the appearance of a massive scalar and 
vector fields describing neutral scalar meson and vector meson, respectively [Ej. The ra- 
diation corrections to the Lagrangian of the massless scalar electrodynamics [H] resembles 
very much, in particular, in their mathematical form at D = 4, the magnetic fluctuation 
corrections to the GL free energy of 4D superconductors [Tj. The interrelationships be- 
tween the superconductivity theory and the gauge theories of elementary particles has 
been comprehensively discussed in Ref. [15]. Note also the interrelationship between the 
GL functional of superconductivity and the CP^^^ confinement model (see Ref. [16]) and 
extensions to non-Abelian theories [17] . 

The gauge theories, mentioned so far, and their extensions have a wide application in the 
description of the Early Universe [T^ [T5] . Another interesting gauge theory is that of 
the nematic-to-smectic A phase transition in liquid crystals. According to the Kobayashi- 
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McMillan-de Gennes theory [1^1 [201 EI] the smectic-74 order is described by two order 
parameters: the nematic director vector and the complex scalar describing the center of 
masses of the long molecules. When a description, quite analogous to that of supercon- 
ductors, is introduced, as suggested by de Gennes [2l|, the director vector is substituted 
with a gauge vector field, which is quite similar to the vector potential A in the GL func- 
tional. Apart from some specific features intended to take into account the liquid crystal 
anisotropy, the effective free energy of the smectic A looks very like the GL free energy 
of superconductors (see, also, Refs. [221 ISSl [21]). Another gauge theory in condensed 
matter physics that has some (but not very close) similarity with the GL free energy 
of superconductors is the Chern-Simons-Ginzburg-Landau (CSGL) effective functional of 
the quantum Hall liquid state, in particular, in the description of the phase transitions 
between the plateaus in the quantum Hall effect [25l |26l [271 [2H]- Let us mention also the 
liquid metallic hydrogen, where the problem of the superconducting-to-superfluid phase 
transition [29] is also related to the topics, discussed in this review. 

1.3 About the investigation of the fluctuation-driven flrst order 
phase transition in superconductors 

Now we shall focus on the fluctuation effects on the properties of the normal-to-superconducting 
phase transition in zero external magnetic field, which is the subject of the present review 
and in particular cases we shall refer to related topics in other natural systems. We shall 
consider in more details the HLM effect in 2D and quasi-2D superconductors with a spe- 
cial emphasis on the problem of the theoretical predictions reliability. Following Ref. [6] 
we shall use two theoretical methods. 

In Sec. II we use the MF like approximation of Ref. [6] where the spatial fluctuations 
of superconducting order parameter ip are neglected, and the magnetic fluctuations are 
taken into consideration. This treatment was justified for well established type I super- 
conductors, where the London penetration length A is much smaller than the coherence 
(correlation) length ^. Due to the neglecting of spatial fluctuations of ip, the results of such 
MF treatment should be valid outside the Ginzburg critical region ^ . A weakly first order 
phase transition occurs as a result of new small |?/^|^-term, which appears in the effective 
free energy of 3D superconductors within the framework of the MF like approximation [6] . 
But this HLM effect was found to be very small and experimentally unobservable even 
for well established type I bulk (3D-) superconductors, such as Al, where the GL number 
K = X/^ = 10-2 < 1. It has been recently shown [301 [3ll [321 [331 [Ml [351 [36] that the 
HLM effect is much more strong in quasi-2D superconductors than in bulk (3D) samples. 
Moreover, as shown in this series of papers, the effect appears by a term of type |'?/'|2ln|'?/'| 
in the effective free energy, and essentially depends on the thickness Lq of quasi-2D su- 
perconducting films. These circumstances provide a real opportunity for an experimental 
verification of the effect - a topic of discussion throughout the present review, and in 
particular, in Sec. II. 
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We must emphasize that MF results are not vahd for thermodynamic states in the 
Ginzburg critical region {ST)g = \Tg — Tco\/Tco around the equilibrium phase transi- 
tion point Tco- In certain classes of high-temperature superconductors the Ginzburg re- 
gion exceeds O.IK whereas in usual low-temperature superconductors it is very narrow, 
{ST)g ~ 10^^^ — 10^^^ K, and the respective critical states are experimentally unaccessi- 
ble (see, e.g., [3]). If the metastability states, which go along with the weakly first order 
phase transition predicted by MF, extend over temperature intervals larger than the size 
of the Ginzburg region, one may conclude that the MF prediction of the HLM effect is 
reliable. In Sec. II we justify the MF analysis reliability for element superconductors as 
Al, W, In. We will show that the condition for enough wideness of metastability regions 
is quite strong and, perhaps, irrelevant in experiments. We will also demonstrate that the 
magnitude of the critical magnetic field is crucial for the observability of the HLM effect 
in thin superconducting films (see also the discussion in Sec. II. 9). 

In Sec. Ill we review some results for the HLM effect obtained with the help of the 
renormalization group (RG) method [U [5]. The latter allows simultaneous treatment 
of both superconducting and magnetic fluctuations in the asymptotically close vicinity 
of the phase transition point. The lack of fixed point of the one-loop RG equations 
for conventional superconductors in zero external magnetic field was interpreted as a 
signal for a fluctuation-driven first order phase transition [6]. This result shows that 
the local gauge magnetic fluctuations are relevant also in the Ginzburg region of strong 
'(/'-fluctuations and under certain circumstances they can change the order of the phase 
transition to a weakly first order, as is outside the critical region. The RG investigations 
of the superconducting phase transition order in zero magnetic field has been recently 
reviewed in Ref. [3^. While the latter review emphasizes the investigation of the HLM 
effect in high orders of the loop expansion [H E], here we lay stress on RG results for the 
effects of anisotropy, quenched disorder and extended symmetries of the Higgs ^-field on 
the phase transition properties. The RG study of unconventional superconductors will 
also be discussed. In the remainder of this paper we shall consider mainly the phase 
transition to the superconducting state but results for other systems and related topics 
will also be mentioned in brief; see Sec. HI. 3 - HI. 5. 



2.1 General GL functional 

The GL free energy [2] of D-dimensional superconductor of volume Vd = {Li...Lci) is 
given in the form 
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In Eq. (1) the first Landau parameter a = ao{T — T^o) is expressed by the critical tem- 
perature Tco = Tc{H = 0) in zero external magnetic field {H = \H\), 6 > is the second 
Landau parameter and e = |e| is the electron charge. The square of the magnetic 
induction B = {H + inM), is given by the vector potential A{x) = {Aj{x), j = 1, D} 
in the form 



B 



2 ^ \ dxi dxj 



here the vector potential A{x) obeys the Coulomb gauge V ■ A{x) = 0. For 3D supercon- 
ductor the relation B = V x A{x) can be used and when B = Bq is uniform along the z- 
axis, the Landau gauge Ao{x) = Bo{—y/2,—x/2,0) can be applied. This representation 
can be generalized for > 2 - dimensional systems, where the magnetic induction Bq is 
second rank tensor |38j : 

Boij = Bo{5ii5j2 — Sj2Sii). (3) 

If we use the notation x = {xi,X2,r), where r is a — 2) - dimensional vector, perpen- 
dicular to the plane {xi, X2), in the 3D case we will have r = (0, 0, z), and 

= 2^jkiBoki = BqSj^ , (4) 

where ejki is the antisymmetric Levi-Civita symbol. The Landau gauge and Eqs. ([3]) - (jlj) 
can be used for uniform B = Bq when 5B - fiuctuations are neglected. In the prevailing 
part of our study we shall apply the general Coulomb gauge of the field which does 

not exclude spatial dependent magnetic fluctuations 5B{x). 

In nonmagnetic superconductors, where the mean value (M) = (M — 5M) of magnetiza- 
tion M is equal to zero in the normal state in zero external magnetic field, the magnetic 
induction in presence of external magnetic field takes the form: 

B = Ho + 6H{x) + 47t6M{x) , (5) 

— * — * 

where Hq is the (uniform) regular part of the external magnetic field and SH is an irregular 
part of H created by uncontrollable effects. We neglect the irregular part 5H and set 
Hq = 0, then B contains only a fluctuation part B = 6B{x) = Att6M{x) that describes 
the diamagnetic variations of M{x) around the zero value (M) = due to fluctuations 
Sijj{x) of the ordering fleld ip{x) above (T > T^o) and below (T < Tcq) the normal-to- 
superconducting transition at Tcq. Note, that the non-fluctuation part Aq = [A{x)—6A{x)] 
corresponds to the regular part Bq = {Hq+{M)) = of i? in nonmagnetic superconductors 
((M) = 0) in zero external magnetic fleld {Hq = 0). Then we can set Aq{x) = and, 
hence, 6A{x) = A{x), so we have an entirely fluctuation vector potential A(x), which 
interacts with the order parameter iIj{x). This interaction can be of type and 
I'j/'pA^ and generates all effects discussed in the paper. 

We accept periodic boundary conditions for the superconductor surface. This means 
to ignore the surface energy including the additional energy due to the magnetic fleld 



6 



penetration in the surface layer of thickness equal to the London penetration depth A(T) = 
-^ol^ol""^^^! ^0 = \T — Tcq\/Tcq] Aq = (mc^fe/Syre^aoT'co)"'^''^ is the "zero-temperature" value 
of A. This approximation is adequate for superconductors of thickness Lq ^ A(T) ^ 
ao, where qq is the lattice constant and Lq = min{Li, i = 1,...,D}. As we suppose 
the external magnetic field to be zero {Hq = 0) or very small in real experiments, the 
requirement Lq ^ A(T) can be ignored and we have the simple condition Lq ^ a^. 

In microscopic models of periodic structures the periodic boundary conditions confine the 
wave vectors ki = {ki = {2'Kni/Li); i = 1,...,D} in the first Brillouin zone [— (vr/ao) < 
ki < (vr/ao)] and the expansion of their values beyond this zone can be made either by 
neglecting the periodicity of the crystal structure or on the basis of the assumption that 
large wave numbers k = \k\ have a negligible contribution to the calculated quantities. 
The last argument is widely accepted in the phase transition theory, where the long- 
wavelength limit {kaQ <^ 1) can be used. In particular, this argument is valid in the 
continuum limit {Vd/clq oo). Therefore, for both crystal and nonperiodic structures 
we can use the cutoff A ~ (7r/ao) and afterwards extend this cutoff to infinity, provided the 
main contributions in the summations over k come from the relatively small wavenumbers 
(A; ^ A). This is in fact a quasimacroscopic description based on the GL functional (1), 
which means that the microscopic phenomena are excluded from our consideration. 

The GL free energy functional takes into account phenomena with characteristic lengths 
.^0 and Ao or larger and A), where A(T) is the London penetration length mentioned 
above and ^(T) = ^oN"^''^ is the coherence length [2]; here ^o = (^^/4mQ;oTco)"'^''^ is the 
zero-temperature coherence length. In low-temperature superconductors ^o and Ao are 
much bigger than the lattice constant ao. Having in mind this argument we will assume 
in our investigation that A <^ (vr/ao). Whether the upper cutoff A is chosen to be either 
A ~ 1/^0 or A ~ 1/Ao is a problem that has to be solved by additional considerations. 
According to arguments presented in Ref. [30] and Sec. II. 6, we will often make the choice 

We will use the Fourier expansion 

M^) = -Lj2Mk)e''-' (6) 

k 

and 

k 

where the Fourier amplitudes Aj(k) obey the relation A*{k) = Aj(—k) and k.A{k) = 0. 
The Fourier amplitude ipi^) is not equal to ip*{—k) because ip{x) is complex function. 
For the same reason ip{0) = tp{k = 0) is complex number. 

The functional (1) is invariant under global U{1) rotations, defined by iplx) —>■ ^/'(x)exp(zQ;), 
where the angle a does not depend on the spatial vector x, and under the local f/(l) gauge 
transformations ip{x) il!{x)exp[ia{x)], A{x) — A{x) + {hc/2e)Va{x). According to the 
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discussion in Sec. I, we have to investigate the spontaneous breaking of these symmetries, 
that is, the ordered phases and the phase transitions in the superconductor. The HLM 
effect, on which we are going to focus the attention, is one of the results of the local U{1) 
symmetry breaking. 

2.2 Notes about the MF like approximation 

While the effect of the superconducting fluctuations Sipix) on the phase transition prop- 
erties is very weak, as in usual superconductors, and is restricted in a negligibly small 
vicinity (|to| ~ 10^^^ ^ 10^"'^^) of temperature Tco, we will assume that 6ip{x) = 0, i.e., 
ip ^ (iplx))] from now on we will denote {ip{x)) by tp. So we apply the mean-field approx- 
imation with respect to the order parameter iplx). Within this approximation we will 
take into account the 5y4(x)-fiuctuations for Bq = 0, i.e., A{x) = 6A{x). Furthermore, 
the y4(x)-fiuctuations can be integrated out from the partition function, defined by: 

Z{^) = J vAe-^^^'^^""^^ , (8) 

where the functional integral j T>A is given by 

oo D 

n n dA,ix)5[^vAix)] . (9) 

The integration is over all possible configurations of the field A{x); the 5-function takes 
into account the Coulomb gauge. 

The partition function Z{iIj) corresponds to an effective free energy jF-p of the dimensional 
system: 

= -kBT\nZ{iP). (10) 

The magnetic fluctuations will be completely taken into account, if only we are able 
to solve exactly the integral ([8]). The exact solution can be done for a uniform order 
parameter ip. The uniform value of ip is different from the mean- field value of ip, because 
the uniform fluctuations of ip{x) always exist, so we should choose one of these two 
possibilities [SUl ES]- The problem of this choice arises after calculating the integral dHD 
at the next stage of consideration when the effective free energy jF-p is analyzed and 
the properties of the superconducting phase {ip > 0) are investigated. The effective free 
energy is a particular case of the effective thermodynamic potential in the phase transition 
theory [H [S] and we must treat the uniform ip in the way prescribed in the field theory 
of phase transitions. It will become obvious from the next discussion that we will use a 
loop-like expansion, which can be exactly summed up to give a logarithmic dependence 
on \ip\. 

Due to the spontaneous symmetry breaking of the global U{1) continuous symmetry of 
the ground state ip ^ 0, the effective free energies discussed in this Section depend on the 
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modulus {ipl of the complex number = j-^le*^ but not on the phase angle 6, which remains 
arbitrary. That is why we will consider the modulus as an "effective order parameter" 
as the angle 6 does not play any role in the phenomena investigated in this Section. 
The quantity remains undetermined up to the stage when we define the equilibrium 
order parameter \ipo\ by the equation of state [dJ^D{4')/d4'] = 0. This equation gives the 
equilibrium value ipo of and the difference 6ipo = (V'o - V') can be treated as the uniform 
(zero dimensional) fluctuation of the field ip{x). The af-dependent fluctuations 6ip{x) have 
been neglected because of the uniformity of ip. The solution ipQ will be stable towards the 
uniform fluctuation 6ip, provided the same solution ipo = \ipo\e''^° corresponds to a stable 
(normal or superconducting) phase; the phase angle 6o remains unspecified. We begin 
our investigation by setting ip uniform but at some stage we will also ignore the uniform 
fluctuation Sip and deal only with the equilibrium value i/jq oi ip. The equilibrium value 
will be calculated after taking into account magnetic fluctuations, so it will be different 
from the usual result \ipo\ = (kl/^)^^^ [2] where both magnetic and superconducting 
fluctuations are ignored. This simplest approximation for the equilibrium value of ip is 
obtained from the GL free energy (1), provided e = and the gradient term is neglected. 
Hereafter we will keep the symbol \ipo\ for the equilibrium order parameter in the more 
general case when the magnetic fluctuations are not neglected and will denote the same 
quantity for e = by ?7 = \ipo{e = 0)| = (|a|/6)^/^. 

The above described approximation neglects the saddle point solutions of GL equations, 
where {ipix)) is x-dependent. Therefore, the vortex state that is stable in type II super- 
conductors cannot be achieved. This is consistent with setting the external magnetic field 
to zero, so the vortex state cannot occur in any type superconductor. These arguments 
can be easily verified with the help of GL equations [2] for zero external magnetic field; the 
only nonzero solution for in this case is given by r] = (|a|/6)^/^ although the magnetic 
fluctuations A{x) = 6A{x) are properly considered. 

In conclusion we can argue that the described method will be convenient for both type I 
and type II superconductors in zero external magnetic field, if the V'-fluctuations have a 
negligibly small effect on phase transition properties Tco = Tc{Hq = 0), where Tc denotes 
the phase transition line for any Hq > 0. For type II superconductors in Hq > 0, two 
lines Tci{Hq) and Tc2{Hq) should be defined, usually given by HdiT) and Hc2{T) |2]. 

2.3 Effective free energy 

When the order parameter ip is uniform the functional (1) is reduced to 



F{,P,A)=F,{,P) + Fa{iP) 



(11) 



with 



Fo(^) = v^D(ai^r+^i^r) 



(12) 
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and 



I ij=l ^ 



dxi dxj 



(13) 



Here p = PolV'P and po = (Svre^/mc^). It is convenient to calculate the partition function 
Z{ip) and the effective free energy J^Bi"^) in the fc-space, where Eqs. and (fT5l) take the 
form 



D fc<A 



and 
Here 

and 



Y[ Y[ dReAj{k)dImAj{k)6 k ■ A{k) 

FaW = Fa{0) + AFaW . 
FA{0) = ^J2^'\A,ik)' , 



(14) 

(15) 
(16) 

(17) 



AFAm=pJ2\Mk) 

note, that we have used the Coulomb gauge k.A{k) = 0. 
Then the partition function ([H]) will be 

Z{iIj) = e-^«(^)/*^^^ZA(V') , (18) 

where 

Za(^) = J VAe-^^^^^^''^^ (19) 

with FaO') given by (fT5|) and the functional integration defined by the rule ( IT^ . With 
the help of Eqs. (ITOl) - (fT9|) the effective free energy J^d{jP) becomes 



^b(^) = Fo(V^)+^/(V') 
where Fq{%Ij) is given by Eq. (fT2|) and 

■Z(^) 



Z(0)J 



(20) 



(21) 



is the '?/^- dependent fluctuation part of T{ip)- In Eq- (EHD the ^'-independent fluctuation 
energy {—ksT In [2^(0)]} has been omitted. This energy should be ascribed to the normal 
state of the superconductor, which by convention is set equal to zero. 



Defining the statistical averages as 



((■■•)) 



(22) 
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we can write Eq. (12 ip in the form 



(23) 



Eq. ( l23l) is a good starting point for the perturbation calculation of Tfi^ip). We expand 
the exponent in Eq. (!23|) and also take into account the effect of the logarithm on the 
infinite series. In result we obtain 

W) = El!(i^(^^AW)- (24) 

where (...)c denotes connected averages Now we have to calculate averages of type 

{A^{ki),Ap{k2)...A^{kn))c. (25) 
Here we will use the Wick theorem and the correlation function of form 

Gf(k, k') = {A,{k)A,{-k')) = 6^,gGf^ik) , (26) 

where 

G0) = {A{k)A,{-k)) = (5,, - hk,) (27) 

and ki = (ki/k). 

The calculation of lowest order terms (/ = 1,2,3) in Eq. ( !24l) with the help of ( l25|) - ( l27l) 
is straightforward. The perturbation terms in ( l24l) are shown by diagrams in Fig. 1. The 
infinite series (1241) can be exactly summed up and the result is the following logarithmic 
function 



1 , M 
k^ 



(28) 



The same result for J^f{ip) can be obtained by direct calculation of the Gaussian functional 
integral ([8]). This is done using the integral representation of 5- function in ([9]) or (HM 
but it introduces an additional functional integration that should be carried out after the 
integration over Aj{x). 

Eqs. ( ITO|l . (!20l) and (!28l) give the effective free energy density 



in the form 
where 

and 



/d(^)=^d(^)/Vb (29) 



/d(^) = /o(V') + A/^(^), (30) 



/o(^) = a|^|2 + ^|^r (31) 
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Figure 1: Diagrammatic representation of the series (24); • represents the p- vertex in (13) 
and (17), and the sohd hnes represent bare correlation functions {\Aj{k)\'^). 

Eqs. (I20l) and (!29|) - (132|) are the basis of our further considerations. We should mention 
that the fluctuation contribution Afoiip) to f{ip) transforms to convergent integral in the 
continuum limit 




where Kjj = 2^^^7r^^/^/r(D/2) for all spatial dimensionalities D > 2. But the terms 
in the expansion of the logarithm in (132!) are power-type divergent with the exception 
of several low-order terms in certain dimensionalities D. Therefore, we will work with 
a finite sum of an infinite series with infinite terms. In our further calculations we will 
keep the cutoff A finite for all relevant terms in Afulip). This is the condition to obtain 
correct results. 

2.4 2D-3D crossover 

The dimensional (2D-3D) crossover has been considered in Refs. [Ml [39] . Here we follow 
Ref. [31], where the effective free energy density f{ip) = fsi^il^) = J^3{4')/V3 of a thin 
superconducting film of thickness Lq and volume V = V^, = (L0L1L2) is derived in a more 
general way, which allows the investigation of the dimensional crossover. Now one should 
perform the integration (33) only with respect to the wave vector components ki, and 
^2 corresponding to the large sizes Lj ^ Lq (j = 1,2). The result for the effective free 
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energy density is 



f{^lj) = a\^\^ + -\i,\' + kBTJ[p{^)] 



where 



is given by the sum 



and q = q= {h, 



J{P) 



S 



1 



+Ao 
fco=— Ao 



(34) 



(35) 



(36) 



In Eqs. the integral J{p) and the sum S{q,p) over the wave vector k = (g, k^) 

are truncated by the upper cutoffs A and Aq. The finite cutoff A is introduced for the 
wave number q and Aq stands for k^. 

As our study is based on the quasimacroscopic GL approach the second cutoff Aq should 
be again related to rather than to the lattice constant Oq, i.e. Aq ~ {1/^q), which means 
that phenomena at distances shorter than are excluded from our consideration. We 
will assume that the lowest possible value of Aq is (vr/^o); as is for A, but we will keep in 
mind that both Aq and A can be extended to infinity, provided the main contribution to 
the integral J{p) and the sum 5* come from the long wavelength limit [q^o ^ 1). 

In a close vicinity of the phase transition point Tco from normal {ip = 0) to Meissner state 
> 0) the parameter p ~ l^/^p is small and the main contribution to the free energy f{ip) 
will be given by the terms in S with small wave vectors /c -C A. This allows an approximate 
but reliable treatment of the 2D-3D crossover by expanding the summation over ko in fl36l) 
to infinity - Aq ~ oo. A variant of the theory when Aq is kept finite (A = Aq = vr/^o) can 
also be developed but the results are too complicated [36]. Performing the summation 
and the integration in Eqs. (l5^ - fl5B]) we obtain J(p) = (A^/27rLo)/(p), where 



lip) 



dy In 



sh (I LpAy/^Ty) 
sh (i LoA^) 



(37) 



The integral (1371) has a logarithmic divergence that corresponds to the infinite contribution 
of magnetic fluctuations to the free energy of normal phase (Tco > 0,ip = 0). Such type 
of divergence is a common property of many phase transition models. In the present 
case, as is in other systems, this divergence is irrelevant, because the divergent term 
does not depend on the order parameter ip and the free energy f{ip) is defined as the 
difference between the total free energies of the superconducting and normal phases: 

fW = {fs-fN). 

Introducing a dimensionless order parameter = {'ipf-ipo), where ipo = {aoTco/bY^"^ is the 
value oi ijj a.t T = 0, we obtain the free energy flM|) in the form 



2toip^ + + 2(1 + to)CIipip' 



(St 
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with I^jiip"^) given by Eq. (jSZ]), 



(l/vr/t)^, A = vr/^o; and 



of state {df /dip = 0) we find two possible phases: ipoo 



From the equation oi siaie yuj /uip = uj we miu lwo p 
superconducting phase {ipo > 0) defined by the equation 



^0 + ^0 + 4^^3^2 ^(/^<^o) 



where 



i^(z) 



^ ^th (I LqAv/^T^) 



(39) 
and the 

(40) 
(41) 



'0 VV + z 

The analysis of the stability condition {d'^f/dip'^ > 0) shows that the normal phase is a 
minimum of f{ip) for to > 0, whereas the superconducting phase is a minimum of f{ip) if 



where 



k{z) 



dy 



coth (i LqA^u + 2;) 



/o 2/ + ^ 

The entropy jump As = (AS/F) 



2sh2 (i LoAv/^+7) 



(42) 



(43) 



[— (i/((y9o)/(iT] per unit volume at the equilibrium 
point of the phase transition 7^ T^o is obtained in the form 



As(T,: 



cO 



1 + 



C/(^, 



<^cO 



(44) 



where y^co = V'o(^c) is the jump of the dimensionless order parameter at T^. 



The second term in As can be neglected. In fact, taking into account the equation 
/[v9o(Tc)] = for the equilibrium phase transition point we obtain that |C/((^o)/v'ol 
approximately equal to Itco+V^co/^l; where ^j^q and the dimensionless shift of the transition 
temperature tco = to{Tc) are expected to be much smaller than unity. The latent heat 
Q = TcAs{Tc) and the jump of the specific heat capacity at Tc, AC = Tc{dAS/dT) can 
be easily calculated with the help of Eq. (jH]). For this purpose we need the function 
ipQ{T), which cannot be obtained analytically from Eq. fHU]) . 

Eqs. (38) and (40) can be analyzed numerically. This relatively simple 2D-3D crossover 
formulae can be used in investigations of specific substances by variation of the thickness 
Lq of the films from Lq ^ A~^ {3D system) to aQ < Lq <^ (quasi-2D system) and 
even to a 2D system for Lq = oq. Then one can vary the effective dimension of the 
system D(,ff{LoA) as a function of LqAq from D = 2 to D = 3 |13]. However, from a 
purely calculational point of view we have found more convenient to consider particular 
dimensions of interest separately and then to compare the results in order to demonstrate 
the relevant differences between the bulk {3D) and thin film properties. This approach is 
applied below. 
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2.5 Effective free energy for particular dimensions 



For purely 2D superconductor consisting of a single atomic layer, we can use Eqs. (129|) - 
(!32|) setting D = 2 and calculate Af2{i') with the help of the rule (l33ll : 



57r 



A2 I "I'i A2 
The first term of this free energy can be expanded in powers of 



(45) 



A/.«0^(|^){.ol.f..ol.f.n(^).ff}. (4a) 

Thus we obtain the result from Ref. HOl. This case is of special interest because of the 
logarithmic term in the Landau expansion for f{ip) but it has no practical application for 
the lack of ordering in purely 2D superconductors. 

For quasi-2D superconductors we assume that (27r/A) > Lq ^ oq, where Lq is the 
thickness of the superconducting film and the more precise choice of the upper cutoff 
A -C (l/cto) for the wave numbers ki is a matter of additional investigation [30]. In order 
to justify this definition of quasi-2D system one can use the 2D-3D crossover description 
presented in Sec. 2.4. The summation over the wave number kg = (27rno/-Lo) in Eq. 
(36) cannot be substituted with an integration because Lq ^ Lj and the dimension Lq 
does not obey the conditions, valid for Lj, (j = 1, 2) [HI |42l |43]. Therefore, for such 3D 
system we must sum over fco and integrate over two other components (fci and ^2) of the 
wave vector k (see Sec. 2.4). This gives an opportunity for a systematic description of 
the 2D-3D crossover as shown in Sec. 2.4. In the limiting case of very small thickness 
the 2D-3D crossover theory (Sec. 2.4) leads to a result, which is obtained more simply 
in an alternative way, namely, by ignoring all terms corresponding to /cq 7^ in the sum 
in Eq. ([32D- This corresponds to the supposition that the quasi-2D film thickness can- 
not exceed 27rA. Assuming this point of view the real physical size of the quasi-2D film 
thickness will depend on the choice of the cutoff A. It is certain at this stage that A > ^o, 
because the (quasi)phenomenological GL theory does not account phenomena for size less 
than ^0- The upper cutoff A of wave numbers can be defined in a more precise way at 
next stages of consideration; see Sec. II. 6 - II. 8). 

For a quasi-2D film we have the expression: 

A/(^) = f A/2(^) , (47) 

where Af2{'ip) is given by Eq. fH5|) . 

For bulk (3D) superconductor we obtain: 



A/3(^) ^ ^ 



1^ 1 + + oPol^rA - isPo \^\ arctan 



3 V ; 3— 3-" I 



(45 
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The Landau expansion in powers of {ipl in this form of fsiip) confirms the respective 
results of Refs. [6], [7], moreover it correctly gives the term of type PqI?/'!^, which has been 
considered small and neglected in these papers. 



For 4D-systems A/d(V') becomes 
SkeT 



647r2 



AVoi^r^ + A* i.i (1 + - m' 1" (1 + ^) 



(49) 



The above expression for Af^ltp) can be also expanded in powers of to show that 
it contains a term of the type |'?/'|^ In (y/p^|?/'|/A), which produces a first order phase 
transition; this case is considered in the scalar electrodynamics [H], as mentioned in 
Sec. I. In our further investigation we will focus our attention on 3D and quasi-2D 
superconductors. 

The free energy density A/d ('?/') can be expanded in powers of IV"! but the Landau ex- 
pansion can be done only in an incomplete way for even spatial dimensions. Thus f^i'ip), 
f^lijj), and f{ijj) - the free energy density corresponding to the quasi-2D films, contain 
logarithmic terms, which should be kept in their original form in the further treatment of 
the function A/d('?/') in the Landau expansion. The analysis has been performed [35] in 
two ways: with and without Landau expansion of A/^i (■?/'). These variants of the theory 
are called "exact" theory (ET) and "Landau" theory (LT), respectively [33j. It has been 
shown [35] that these two ways of investigation give the same results in all cases except 
for quasi-2D films with relatively small thicknesses {Lq <^ C,o)- It seems important to es- 
tablish the differences between two variants of the theory because the HLM effect is very 
small and any incorrectness in the theoretical analysis may be a cause for an incorrect 
result. By same arguments one can investigate the effect of the factor T in A/d ('?/') on the 
thermodynamics of quasi-2D films. This factor can be represented as T = Tco(l + to) and 
one may expect that the usual approximation T ^ Tco, which is well justified in the Lan- 
dau theory of phase transitions [2],|1], may be apphed. This way of approximation can be 
made by neglecting terms in the thermodynamic quantities smaller than the leading ones. 
On the other hand practical calculations lead to the conclusion that this approximation 
cannot be made without a preliminary examination because for some quasi-2D films it 
produces a substantial error of about 10% [SS]- LT, in which the factor T is substituted 
by TcQ is referred to as "simplified Landau expansion" - SLT. All three variants of the 
theory, ET, LT and SLT have been investigated in Ref. [35] . 



2.6 Limitations of the theory 

The general result (l29l) - ( l32l) for the effective free energy f{ip) has the same domain of va- 
lidity [2] as the GL free energy functional in zero external magnetic field. When we neglect 
a sub-nano interval of temperatures near the phase transition point we can use Eq. (1), 
provided |to| = \T — Tco\/Tco < 1, or in the particular case of type I superconductors, 
\to\ < 1^"^ i2j- Note, that the latter inequality does not appear in the general GL approach. 
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It comes as a condition for the consistency of this approach with the microscopic BCS 
theory for type I superconductors [2]. 

Taking the continuum hmit we have to assume that all dimensions of the body, including 
the thickness Lq, are much larger than the characteristic lengths ^ and A. The exception of 
this rule is when we consider thin films. Especially for thin films of type I superconductors, 
where ((27r/A) > Lq ^ ao), we should have in mind that ^(T) > A(T), so the inequalities 
,^ > A > ^0 > -^0 hold true in the domain of validity of the GL theory |to| < < 1- 
In Ref. [30] a comprehensive choice of the cutoff A has been made, namely, A = (the 
problem for the choice of cutoff A is discussed also in Sec. II. 7 - II. 8). The respective 
conditions for quasi-2D films of type II superconductors are much weaker and are reduced 
to the usual requirements: k > l/-\/2, |to| < 1 and (27r/A) > Lq ^ uq. 

If we do a Landau expansion of foii')-, in powers of the condition p <^ should 
be satisfied. In order to evaluate this condition we substitute l^/'p in p = polV^P with 
rf = \a\/b, which corresponds to e = 0. As A^(T) = 1/p, the condition for the validity of 
the Landau expansion becomes [AA(T)]^ ^ 1, i. e., (AAq)^ ^ |to|- Choosing the general 
form of A^ = (vrr/^o), where r describes the deviation of A,- from Ai = A = (vr/^o), we 
obtain (vrr/t)^ ^ \to\; n = (Aq/^o) is the GL parameter. 

Thus we can conclude that in type II superconductors, where k = (Aq/^o) > l/v^5 ^^e 
condition (p/A^) -C 1 is satisfied very well for values of the cutoff in the interval between 
A = (vr/^o) and A = (tt/Aq), i.e., for 1 < r < (I/k). For type I superconductors, 
where k < 1/V2 the cutoff value A ~ (l/^o) leads to the BCS condition (|to| < /t^) for 
the validity of the GL approach. Substantially larger cutoffs (A ^ tt/^o), for example, 
A ~ (1/Ao) for type I superconductors with k <ti I lead to contradiction between this 
BCS condition and the requirement p <^ A^. 

In our calculations we often use another parameter p^ = (l/vrr/t)^ and, in particular, 
p = Pi = (I/ttk)^ and in terms of p the condition for the validity of expansion of foi'^P) 
becomes p|to| <^ 1, or, more generally, Prl^ol 1- Choosing r = l/vr we obtain the 
BCS criterion for the validity of the GL free energy for type I superconductors [2j. The 
choice T = (^o/ttAo) corresponds to the cutoff A,- = I/Aq. As we will see in Sec. II. 7 
the thermodynamics near the phase transition point in 3D systems has no substantial 
dependence on the value of the cutoff A^ but it should be chosen in a way that is consistent 
with the mean-field-like approximation. The cutoff problem for quasi-2D films has been 
investigated in Ref. [35]. It has been shown that the choice A ~ tt/^q is consistent in the 
framework of GL theory. 

Alternatively, the inequality (p/A^) ^ 1 may be investigated with the help of the reduced 
order parameter defined hy = \'ip\/r]o, where r]o = r](T = 0) = {aoTco/by^"^ is the so- 
called zero-temperature value of the order parameter for the GL free energy fo{ip), given 
by Eq. fl3T]) . The reduced order parameter (p will be equal to |to| for to < 0, if only the 
magnetic fluctuations are ignored, i.e., when = rj. Using the notation (p, we obtain 
the condition {p/ h?) ^ 1 in the form PrV'^ <^ 1- This condition seems to be more precise 
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because it takes into account the effect of magnetic fluctuations on the order parameter 



2.7 Bulk superconductors 
2.7.1 Thermodynamics 

The effective free energy fsi^ip) of bulk (3D-) superconductors is given by Eqs. (l29l) - (!3Ti) 
and (148|) . The analytical treatment of this free energy can be done by Landau expansion 
in small {y/polipl / A) . Up to order \ip\^ we obtain 

Mi^) ^ aslV^P + - qM' + flV^r , (50) 



where 



ksTApo 
27r2 

27r2A 
Qn 



03 = aH — — , (51) 

(52) 

.3 = ^ , (53) 



and 



- "6^ • ^^^^ 
The cutoff A in Eqs. (15T]) - (15^ is not specified and can be written in the form A,- = (vrr/^o) 
as suggested in Sec. 2.6. 

It can be shown by both analytical and numerical calculations [31] that |^|^-term has no 
substantial effect on the thermodynamics, described by the free energy flH?]) . That is why 
we ignore this term. The possible phases \ipo\ are found as a solution of the equation of 
state: 

[dfW/dm = . (55) 

There always exists a normal phase {ipol = 0, which is a minimum of /s ('?/') for 03 > 0. 
The possible superconducting phases are given by 




160363 
9g| 



|V>o|± = ^ I 1 ± Jl - 1 > 0. (56) 



Having in mind the existence and stability conditions of |?/'o|±-phases jSj, we obtain that 
the l^/'ol+'Phase exists for (160363) < 9g| and this region of existence always corresponds 
to a minimum of fsiip)- The |'?/'o|--phase exists for < 03 < 9g|/1663 and this region of 
existence always corresponds to a maximum of fs^ip), i.e., this phase is absolutely unstable. 
For 03 = 0, iV'ol- = and hence, it coincides with the normal phase. For 9g| = 160363 
we have |^o|+ = \i'o\- = ^Qs/^h and f3{\i'o\+) = fsHi'ol^) = 27g|/5126|. Furthermore 
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/sdV'ol-) > for all allowed values of lipol- > 0, whereas /3(|V'o|+) < for 03 < gl/263, 
and fsdipol^) > for (gl/263) < 03 < 9^3/1663. The equilibrium temperature Teq of 
the first order phase transition is defined by the equation /(|?/^o|+) = 0, which gives the 
following result: 

263(Teq)a3(Teq) = ql{%^) . (57) 

These results are confirmed by numerical calculations of the effective free energy (ISOll [31] . 

The equilibrium entropy jump is AS* = VAs and As = —{df^dipD/dT) can be calculated 
with the help of Eq. ( l50l) and the equation of state ( !55l) : 

As = -|^or$(|^o|), (58) 
where ^(I'j/'oI) is the following function: 

The specific heat capacity per unit volume AC = T{dAs/dT) is obtained from (!58l) : 

AC=-f^«t(W). (60) 



0\ 

Tco J dto 

The quantities As{T) and AC(T) can be evaluated at the equilibrium phase transition 
point Teq, which is found from Eq. fl57j) : 

Teo 27r2ao & + (pgfcB/27r2A) T^o V "o 7 ' ^ ^ 

provided |ATc| = \Tco — Teq| <C Tcq. Further we will see that the condition |ATc| <^ T^q 
is valid in real substances. The second term in r.h.s. of Eq. ( l6TI) is a typical negative 
fluctuation contribution whereas the positive third term in r.h.s. of the same equality is 
typical of first-order transitions [1]. 

To obtain the jumps As and AC at Tgq we have to put the solution |'?/'o|+ found from 
Eq. ([56D in Eqs. ([58]) - ([60]). The result will be: 



and 

4ao / ^ qlb 



AC = ^^.„T,-^j, (63) 

where 63c and q^c are the parameters 63 and 53 at T = Teq. As |ATc| = \Tcq — Teq| ^ T^o 
we can set Teq ~ Tco in r.h.s. of Eqs. (l62l) and (!63l) and obtain q^c = qz{T = Teq) ~ g3(Tco) 
and 63c ^ 63 (^co)- 
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The latent heat Q = — TgqAs of the first order phase transition at Tgq can be calculated 
from Eq. fl62|) . If we neglect the second and third terms in the r.h.s. of Eq. (62) and the 
second term in the r.h.s. of Eq. (63) we shall obtain a result for the ratio 

(AT)eq = ^ . (64) 

which is similar to that presented in Ref . [6] . Here we should mention that Eq. fl63|) gives 
the jump AC at the equilibrium phase transition point of the first order phase transition, 
described by \ipf term while the specific heat jump considered in Ref. |6j is equal to 
the specific heat jump at the standard second order transition AC = {alTco/b)^ and is 
four times smaller. Therefore, we obtain (AT)cq four times smaller than the respective 
value in Ref. [6]. 



2.7.2 Results for Al 



In order to do the numerical estimates we represent the Landau parameters ao and b with 
the help of the zero-temperature coherence length and the zero-temperature critical 
magnetic field Hcq. The connection between them is given by formulae of the standard 
GL superconductivity theory \T\: C,q = (/i^/4mao^co) and H'^q = (inalT^Q/b). The 
expression for the zero-temperature penetration depth Aq = (^c/2v^eifco^o) is obtained 
from the above relation and Aq = (6/ao^coPo)^''^- We will use the following experimental 
values of T^o, Hco and for Al given in Table 1. The experimental values for T^o, H^o 
and ^0 vary about 10-15% depending on the method of measurement and the geometry 
of the samples (bulk material or films) [S] but such deviations do not affect the results 
of our numerical investigations. 



Table 1. Values of T^o, Hcp, ^p, k, and [V^qI for W, Al, In [B]. 



substance 


Tco (K) 


Hco (Oe) 


Co (^m) 




IV'ol X 10-" 


W 


0.015 


1.15 


37 


0.001 


0.69 


Al 


1.19 


99.00 


1.16 


0.010 


2.55 


In 


3.40 


281.5 


0.44 


0.145 


2.0 



The evaluation of the parameters 03 and 63 for Al gives: 

03 = (aoTco) [to + 0.972 x 10-^(1 + to)T] , (65) 

and 

- = 1 + -^. (66) 

Setting T = 1 corresponds to the cutoff Ai = (7r/,^o) (Sec 2.6). For r = (1/k)^; = 10^, which 
corresponds to the much higher cutoff A = (tt/Aq) we have 63 « b. i.e., the Pq -term in 63, 
given by Eq. (j52p . can be neglected. However, as we can see from Eq. (j66p . for r = 1 the 
same pQ-correction in the parameter 63 is of order 0.16 and cannot be automatically ignored 
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in all calculations, in contrast to the supposition in Refs. OE]. However, the more important 
fluctuation contribution in 3D superconductors comes from the r— term in Eq. (j65p for the 
parameter 03. This term is of order 10~^ for r ~ 1 and this is consistent with the condition 
1^0 1 < ~ 10~^ but for r ~ 10^, i.e., for A ~ (vt/Aq) ~ lO^/xm, the same r— term is of order 
10^, which exceeds the temperature interval (Tco ± 10~^) for the validity of BCS condition for 
Al (Sec. 2.6). 

These results demonstrate that for our theory to be consistent, we must choose the cutoff 
A,- = (vrr/^o)) where r is not a large number (r — > 1 10). To be more concrete we set 
A = Al = (tt/^o) as suggested in Ref. [30] . 

The temperature shift teq = io(T'eq) for bulk Al can be estimated with the help of Eq. (f6T]l . We 
obtain that this shift is negative and very small: teq ~ — 10~^. Note, that the second term in the 
r.h.s. of Eq. ()6ip is of order 10~^, provided A ~ (l/'^o) whereas the third term in the r.h.s. of 
the same equality is of order 10~^. Once again the change of the cutoff A to values much higher 
than (vr/^o) will take the system outside the temperature interval, where the BCS condition for 
Al is valid. Let us note, that in Ref. [31] the parameter t corresponds to our present notation 
to- But the numerical calculation of the free energy function f^{ip) in Ref. [31] was made for 
the SLT variant of the theory and the shifted parameter (to + 0.972 x 10~^) was incorrectly 
identified with t and this lead to the wrong conclusion for its positiveness at the equilibrium 
phase transition point Tgq. As a matter of fact, the shifted parameter (to + 0.972 x 10"*^) is 
positive at Teq but teq = to(req) is negative, as firstly noted in Ref. [35] . 

Having in mind these remarks, when we evaluate As and AC for bulk Al we can use simplified 
versions of (I62p and ()63p which means to consider only the first terms in the r.h.s and to take 
93c ~ '73 and 63c ~ 6 at Tcq- In this way we obtain 

Q = -T,oAs = 0.8 X 10^2 r ^^g 1 ^ (67) 

LK . cm'^ J 



and 



AC = 2.62 X 10' 



erg 
Lcm'^ 



The results are consistent with an evaluation of AC for Al as a jump (AC = agTco/6) at the 
second order superconducting transition point [6] that, as we have mentioned above, is four 
times smaller than the jump AC given by Eq. (I68p . 



A complete numerical evaluation of the function f'i{ip) and the jump of the order parameter at 
Teq for bulk Al was presented for the first time in Ref. [31] . The results there confirm that the 
order parameter jump and Q for bulk type I superconductors are very small and can hardly be 
observed in experiments. 

We will finish the presentation of bulk Al with a discussion of the ratio (j64p . It can be also 
written in the form _^ 

(AT), - ^^"^ (^bTco\ ( 



and it differs by a factor 1/4 from the respective result in Ref. [6]. As mentioned in Sec. 2.7.1, 
this difference is due to the fact that we take AC as the jump at the first order transition 
temperature Teq, as given by the first (leading) term in Eq. (63), while in the paper [6j the 
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authors use a hypothetic specific heat jump (AC) at the standard second order phase transition 
point. Prom Eq. (j69p we obtain 



(Ar)eq = 6.7 X W-''{T''H'^,eo), (70) 

and multiplying the number coefficient in the above expression by 4 we can obtain Eq. (10) 
presented in Ref. [6]. 



2.8 Quasi-2D films 

Following Refs. [32l[33] we can present the free energy density /(v^) = {F{Tp)/ L1L2) in the form 

f{ip) = ^ [2to^^ + ip'' + C{1 + to)T{f,^^)] , (71) 

where 

r(y) = (l+y)ln(l + y)-2/lny, (72) 

To obtain Eqs. ([7T]) - ([72|) we have set A = (tt/^q) and introduced the notation </? = \tp\/r]Q; 770 
is defined in Sec. 2.5. Some of the properties of free energy (j7ip were analyzed in Ref. [32j for 
Al films and in Ref. [33] for films of Tungsten ( W) , Indium (In) , and Aluminium ( Al) . Here we 
will briefly discuss the main results. 

The equilibrium order parameter lpq > corresponding to the Meissner phase can be easily 
obtained from the equation df{ip)/dip = and Eq. (fTTI) : 



to + ^t>+ \ In ( 1 + — ) = . (73) 



The logarithmic divergence in Eq. ([73]) has no chance to occur because 990 is always positive and 
does not tend to zero. 

The largest terms in the entropy jump bs and the specific heat jump bC at the equilibrium first 
order phase transition point Tcq are given by [32l [33] 

where (/Jgq = = Teq), and 

= . (75) 

The latent heat of the phase transition [4j is given by Q = —Taabs and Eqs. ([73]) - ([74]) . Since the 
temperatures Teq and Tco have very close values, the difference between the values of Q, bs^ and 
bC at Tco and Tgq, respectively, can also be ignored, for example, |(^C(Teq) — (5C(Tco)|/5C(rco) <C 
1 and we can use either bCiTco) or (5C(Tcq) 



The equations ([7ip - ([73p corresponding to quasi-2D films are quite different from the respective 
equations for bulk (3D-) superconductors but it is easily seen that the relatively large values 
of the order parameter jump in thin films again correspond to relatively small values of the 
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GL parameter k. That is why we consider element superconductors with small values of k and 
study the effect of this parameter, the critical magnetic field Hco and the film thickness Lq on 
the properties of the fluctuation-induced first order phase transition. 

We use experimental data for Tco, -ffcOi Co and k for W, Al, and In published in Ref. ,44J (see 
Table 1). In some cases the GL parameter k can be calculated with the help of the relation 
K, = (Ao/Co) and the available data for and Aq. In other cases it is more convenient to use the 
following representation of the zero-temperature penetration depth: 

2V2eiicoCo 

The value of IV'ol in Table 1 is found from 

/ m \ 1/2 , ^ 

The order parameter dependence on the reduced temperature difference to is shown in Fig. 2 
for Al films of different thicknesses. It is readily seen that the behavior of the function (/9o(to) 
corresponds to a well established phase transition of first order. The vertical dashed lines in 
Fig. 2 indicate the respective values of teq = io(^eq)) at which the equilibrium phase transition 
occurs as well as the equilibrium jump (/9o(T'eq) = V'eq for different thicknesses of the film. The 
parts of the (/7o(io)-curves which extend up to to > *eq describe the metastable (overheated) 
Meissner states, which can appear under certain experimental circumstances (see in Fig. 2 the 
parts of the curves on the r.h.s. of the dashed lines). The value of (p^c^ and the metastable region 
decrease with the increase of the film thickness, which shows that the first order of the phase 
transition is better pronounced in thinner films and this confirms a conclusion in Ref. [32j . 

These results are confirmed by the behavior of the free energy as a function of to. We used 
Eqs. (f7T]) - ([73]) for the calculation of the equilibrium free energy /[v?o(*o)]- The free energy for 
Al films with different thicknesses is shown in Fig. 3. The equilibrium points Teq of the phase 
transition correspond to the intersection of the /((/3o)-curves with the to-axis. It is obvious from 
Fig. 3 that the temperature domain of overheated Meissner states decreases with the increase 
of the thickness Lq. 



Table 2. Values of teq, (peq, and Q (erg/cm^) for films of W, Al, and In with different 
thicknesses Lq (/im) |33]. 
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0.20 
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1 


-0.00029 
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-0.00013 


0.006 


0.23 


-0.00032 


0.016 


2.7 X 10" 


-5 


2 


-0.00017 


0.009 


0.06 


-0.00008 


0.004 


0.10 


-0.00021 


0.013 


1.8 X 10" 


-5 



The shape of the equilibrium order parameter v?o(''^o) in a broad vicinity of the equilibrium 
phase transition for thin films [Lq = 0.05/im) of W, Al, and In was found from Eq. (1731) . 
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Figure 2: Order parameter profile f{to) of Al films of different thicknesses: Lq = 0.05 {im 
("-|-"-line), Lq = 0.1 /im (o), and Lq = 0.3 /im (■ 
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Figure 3: The free energy /(to) for Al films of thickness: Lq 
Lq = 0.1 /im (o), Lq = 0.3 fim (■) [33] . 
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Figure 4: Order parameter profile ^plto) of films of thickness Lq = 0.05 fim: W ("-|-"-line), 
Al (o), and In (•) [33]. 

The result is shown in Fig. 4. The vertical dashed lines in Fig. 4 again indicate the 
respective values of tgq = toiTeq), at which the equilibrium phase transition occurs as well 
as the equilibrium jump (y9o(T'eq) = V^eq in the different superconductors. 

The order parameter jump at the phase transition point of In (the In curve is marked 
by points in Fig. 4) is relatively smaller than for W, and Al, where the GL parameter 
has much lower values. The same is valid for the metastability domains; see the parts of 
the curves in Fig. 4 on the left of the vertical dashed lines. It is obvious from Fig. 4 and 
Table 2 that the equilibrium jump of the reduced order parameter ip^q of W has a slightly 
smaller value than that of Al although the GL number k for W has a ten times lower value 
compared with k of Al. Note, that in Fig. 4 we show the jump of (peq, but the important 
quantity is {ipleq = I'ipol^eq- Using the data for Lq = 0.05/im from Tables 1 and 2 we find 
for l^leq the following values: 0.1 x 10^^ for Al, 0.05 x 10^^ for In, and 0.02 x 10^^ for W. 
This result shows that the value of the critical filed Hcq is also important and should be 
taken into account together with the smallness of GL number when the maximal values 
of the order parameter jump are looked for. Thus the value of the order parameter jump 
at the fluctuation- induced phase transition is maximal, provided small values of the GL 
parameter k are combined with relatively large values of the critical field Hcq. In our case 
Al has the optimal values of these two parameters. 

The importance of the zero-temperature critical magnetic field Hco for the enhancement 
of the jumps of the certain thermodynamic quantities at the equilibrium phase transition 
point Tgq becomes obvious from Eqs. (171|) . (ITS]) and (|77|) . Eq. (177|) shows that the order 
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parameter jump lijleq = IV'olv'eq is large for large values of Hco and ^o- This is consistent 
with the requirement for relatively small values of the GL parameter k as shown by 
Eq. fl76l) . Therefore, the unmeasurable ratio Q/6C discussed in Ref. [6] does not depend 
on the value of the critical field Hco but the quantities Q and SC themselves as well as the 
order parameter jump \ip\eq depend essentially on Hco- The values of the reduced order 
parameter jump ip^q for films of Al, In and W of the same thickness have the same order 
of magnitude while the respective order parameter jump {ipleq is one order of magnitude 
higher for Al than for W, as shown above. The effect of the critical magnetic field iJco on 
the latent heat Q is, however, much stronger and, as is seen from Table 2, the latent heat 
Q in W films is very small and can be neglected while in Al and In films it reaches values, 
which could be measured in suitable experiments. This is so because the latent heat is 
proportional to the difference [H^q/Sti ~ &|'?/'oH between the energies of the ground state 
(superconducting phase at T = 0) and the normal state. It is consistent with the fact that 
the fluctuation contribution to the free energy, i.e., the C— term in the r.h.s. of Eq. (17T1) is 
generated by the term of type I'^/'P / d^xA^lx) in the GL free energy. At T = this free 
energy term is also proportional to the mentioned difference between the free energies of 
the ground and normal states. 

The shift of the phase transition temperature tgq = K^cq — ^co)|/^cO) the reduced value ip^q 
of the equilibrium order parameter jump |V^|eq, and the latent heat Q of the equilibrium 
transition are given for films of different thicknesses and substances in Table 2. This table 
shows that the shift of the phase transition temperature is very small and can be neglected 
in all calculations and experiments based on them. The values for yjgq for different Lq 
and those for \^po\ given in Table 1 confirm the conclusion, which we have made for films 
of Al, In, and W with Lq = 0.05/im. The latent heat Q has maximal values for In, where 
the critical field is the highest for the considered materials. 

2.9 Final remarks 

Our analysis shows that the MP studies of the HLM effect have a well defined domain of 
validity for both 3D- and quasi-2D superconductors. Our conclusion is that the MP theory 
of the magnetic fluctuations in superconductors and, in particular, the MP prediction for 
the fluctuation driven weakly first order phase transition in zero external magnetic field 
in bulk superconductors [B] and quasi-2D superconducting films [301 E21 E3] is reliable and 
can be tested by experiments. While the HLM effect in bulk systems is unobservingly 
small, in quasi-two dimensional superconductors this effect is much stronger and might 
be observed with available experimental techniques. 

Our consideration of quasi-2D superconductors is highly nontrivial in view of the relevance 
of the dependence of the effective Landau parameters on the thickness of the films, Lq. 
We have justified this dependence by simple heuristic arguments and by a reliable con- 
sideration of the 2D-3D crossover. In contrast to initial expectations [6] that films made 
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of superconductors with extremely small GL parameter k such as Al and, in particular, 
W will be the best candidates for an experimental search of the HLM effect, our careful 
analysis (firstly published in Ref. [33] definitely gives somewhat different answer. The Al 
films still remain a good candidate for transport experiments, through which the jump 
of the order parameter at the phase transition point could be measured but surprisingly 
the W films turn out inconvenient for the same purpose due to their very low critical 
field Hf.Q. The importance of the critical magnetic field i^co for fhe clearly manifested 
first-order phase transition has been established and discussed, too. Although In has ten 
times higher GL number k than Al, the In films can be used on equal footing with the 
Al films in experiments intended to prove the order parameter jump. Here the choice of 
one of these materials may depend on other features of experimental convenience. As far 
as caloric experiments are concerned, the In films seem to be the best candidate due to 
their high latent heat. 

As shown in this review, experiments intended to search the HLM effect can be performed 
by both type I and type II superconductors. In experiments the external magnetic field 
cannot be completely eliminated. Then vortex states may occur for H = \H\ > Q below 
Tc = Tc{H) < Tco in type II superconducting films and this will obscure the HLM effect. 
Note, that in both type I and type II superconductors the external magnetic field H 
generates additional entropy jump at the phase transition point Tc{H) and this effect 
can hardly be separated from the entropy jump ( TMl) caused by the magnetic fluctuations 
in the close vicinity of T^q. Therefore, in experiments intended to the search of HLM 
effect the external magnetic field should be minimized as much as possible. For quasi- 
2D superconductors, where the HLM effect is relatively strong and the latent heat can 
exceed several ergs, one should ensure external fields less than 10 Oe, or, in more reliable 
experiments, less than 1 Oe [35j. 

The temperature range around Tco, where the HLM effect is significant has been estimated 
in Ref. [6] to a few microdegrees, and later, this estimate has been confirmed in a more 
accurate way [15]. According to our point of view, the possibilities for the observation 
of the effect are related mainly with its magnitude, because the range of temperatures 
where it occurs always exists, even in the critical region of strong T/^-fiuctuations, as the 
RG studies indicate (Sec. III). The review of the results for Al, In, and W in Sec. II shows 
that the metastability domains (of overheating and overcooling) are much larger than the 
Ginzburg critical region. This result confirms the reliability of the MF treatment. As 
one may see from Figs. 2-4 the metastability temperature interval is relatively larger for 
smaller values of the GL number k, but in order to ensure large values of some measurable 
thermodynamic quantities as, for example, the latent heat Q and the specific heat capacity 
C we must choose a material with large critical field He- The experimental verification 
of the order parameter jump can be made by transport experiments. As we see from 
Figs. 2-4 and Table 2, this quantity has maximal values for W, where the parameter k, is 
very small. Having in mind also the large metastability regions in this material, one may 
conclude that W is a good candidate for a testing the HLM effect by transport experiments 
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(measurements of the superconducting currents), provided some specific disadvantages of 
this material (quite low T^o, etc.) from the experimental point of view do not contradict 
to this our conclusion. Another effect, which is relevant to the present discussion, is the 
known variation of the GL parameter k with the variation of the thickness Lq (see, e.g., 
Ref. ^6]. The parameter k < l/\/2 of a type-I bulk superconductor may change up to 
values corresponding to a type-II superconductor with the decrease of the thickness Lq 
below lO^^m. 

In Ref. [17] a very interesting interrelationship between the first order phase transition and 
the vortex fluctuations in type-I 2D superconductors has been analyzed by both analytical 
and numerical methods. A recent investigation [18] of HLM effect in superconducting 
films has been proven wrong in our comment [19]; see, also, the Corrigendum [50] to the 
work [38] . In further papers the same authors and co-workers confirm [STl [52] our results, 
published in Ref. [301 SB 1321 ESl [331 ES] ) and discussed in this Section. Another paper 
has been intended to the treatment of gauge effects in superconductors with the help of a 
Gaussian effective functional [53], which gives results quite near to mean-field ones, and 
is wrong within the present discussion. 

3 RG STUDIES AND OTHER TOPICS 

3.1 Notes about RG 

The RG methods are intended to investigation the strong fluctuation interactions in the 
critical region of second order phase transitions and multi-critical phenomena [31 [5]. In 
studies of complex systems with several orderings or with influence of additional gauge 
fields, as in Eq. (1), usually the fc— space RG methods are used because of their wider 
applicability. The methods differ from each other in specific technical details but all of 
them lead to the same predictions about the phase transition properties. The common 
feature is that the mentioned methods are based on the so-called loop expansion [31 [5l [12] . 
The RG equations are, in fact, infinite asymptotic series, which are truncated at one-, 
two-, or, in rare cases, higher orders in the loop expansion. The total infinite series 
can be summed only in case of trivial (usually exactly soluble) models and in some very 
exceptional cases as, for example, the RG series for interacting real bosons in the quantum 
limit T ^ 154|. 

For the simple 0^— theory the RG series can be derived and analyzed to high orders in the 
loop expansion whereas for more complex models this can be practically done within the 
one- and two-loop approximations. As the RG series are asymptotic, normally, the one- 
and two- loop orders give all important features of the specific system of interest: (i) the 
presence of stable fixed points (FP) of the RG equations and the associated with them 
types of critical behavior, and (ii) the lack of stable FPs and related conclusions about 
the lack of standard critical or multi-critical behavior. The higher orders of the theory 
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are more relevant to investigations of the asymptotic properties of the loop series than for 
obtaining of qualitatively new critical behavior, or, other new qualitative characteristics of 
the system. The latter are reliably obtained within the one- and two-loop approximations. 

3.2 The order of the phase transition 

In Ref. [6j| the simultaneous effect of ip- and A-fiuctuations in the fluctuation Hamiltonian 
(1) has been investigated by one-loop RG and e = (4 — (i)-expansion [1] for the general 
case of n/2-component complex vector field ip = {ipa]C( = 1,...,?t,/2}; note, that this field 
is equivalent to a n-component real vector field. It has been shown [6] that a stable FP 
exists below four dimensions = (4 — e) only for symmetry indices n > 365.9, which are 
far above the real symmetry index n = 2 for usual superconductors and numbers n = 4, 6 
corresponding to superconductors with certain unconventional Cooper pairings For 
n < 365.9 real FP does not exist at all. Besides, the 1/n-expansion has been used for 
the calculation [6j of the critical exponents in the so-called large-n limit (alias "Hartree 
limit") [H[5], where [n > 365.9) the phase transition is definitely of second order. 

According to the usual interpretation of RG results, the lack of stable FP is an indication 
for a lack of standard second order phase transition. But usually, additional (non-RG) 
arguments are needed to determine the actual order of the phase transition. In the 
present case, the lack of stable FP could be a result of the same mechanism that produces 
a fluctuation- driven weakly first order phase transition in MF approximation (Sec. 2). As 
RG takes into account both fluctuations of ip and A, one can conclude that the result 
for the weakly-first order phase transition in zero external magnetic field is valid for both 
type I and II superconductors. Of course, the size of the effect (the size of jumps of 
energy, latent heat, order parameter) will depend on the specific substance. This path of 
investigations has been further developed in the paper [7], and several problems opened 
in the short Ref. 0, have been solved. 

In Ref. [10], where both 2D and 3D superconductors are considered, the HLM effect has 
been confirmed in one-loop order (annealed disorder [10]) as well as the availability of 
second order phase transition for largen has been proven in another variant of the Hartree 
limit {n oo). 

Next, the effect has been confirmed in a RG investigation of the equation of state below 
the phase transition point T^q ^^\57\. These investigations show that the effect may occur 
under quite restricted conditions for the vertex parameters e and b of the functional (1), in 
particular, the effective charge |e*| exceeds some value. A strong restriction on the HLM 
effect, but in an opposite direction - a requirement for a sufficiently small effective electric 
charge |e*|, has been found also by a Monte Carlo (MC) study [58] of a lattice version of 
the model (1). Another MC study concludes that at fixed effective charge (|e^| = 5) 
the HLM effect strongly depends on the GL parameter k = X/^: it is well established for 
K <^ 1, becomes very weak for k ~ 1 and vanishes for large n. On the basis of this picture 



29 



a proposal is made (SH] for the existence of a tricritical point at some "tricritical" value of 
K. As K is related with the charge e, this proposal seems to be in a conformity with other 
investigations. Note, that the first prediction for change of the phase transition order with 
the variation of k and for the possibility of "tricritical point" to appear at some value 
of K, has been made in Ref. [60] on the basis of duality arguments [61] and analytical 
calculations; see further development of this approach in Refs. [621 [63l [BU [65] . A further 
MC study [66] of type II superconductors indicates that these systems undergo a second 
order phase transition of universality class 3D rather than a weakly first order phase 
transition. The two-loop RG studies [671 168] show that the account of higher orders of 
the loop expansion allows a lowering of the critical value ric = 365.9 but the RG equations 
still have no real FP at D = 3 and n = 2 (conventional 3D superconductors), except for 
a treatment by Pade analysis [68j (see, also, the review [37]). In 2 + e dimensions, the 
nonlinear a model exhibits a second-order phase transition for all values of n > [69] . 
Owing to this result one may suppose that the critical value ric vanishes at some dimension 
2 < c? < 4, as mentioned in Ref. [70]. RG studies [711 E2] at fixed dimension d have also 
been carried out with the aim to determine the phase transition order (see, also. Ref. [70] 
for a comment). 

3.3 Disorder effects 

Here we will discuss mainly disorder described by Gaussian distributions [1] . The effects of 
annealed and quenched disorder in classical versions of Abelian-Higgs models, equivalent 
to the GL model (1) have been investigated by Hertz [731 [71] in the context of the theory 
of spin glasses in case of Dzyaloshinskii-Moriya interaction. The specific feature of this 
approach is that the disorder is associated with the vector gauge field A{x) and can be 
used to describe superconductors in random, uncorrelated, or, which is the same, short- 
range correlated magnetic fields rather than usual ones. In high energy physics this 
approach makes contact with gauge fields where a Higgs field is coupled to a random color 
field. It is natural to expect, as has been proven [10l[73], that the annealed gauge model 
will bring a fluctuation-driven flrst order phase transition at usual symmetry indices and 
a second order phase transition in the Hartree limit (n — > oo). In Ref. [40] the phase 
transition in case of quenched impurities has been predicted to be of second order in 
a calculation to one-loop order, in contradiction to the conclusion for a lack of stable 
nontrivial FP of the RG equations within the same one-loop order [73] . The origin of this 
discrepancy is clear from the argument that the lack of stable FP may result without any 
change of the Hamiltonian structure, whereas the conclusion for a second order of the 
transition in Ref. [ID] has been made, perhaps, incorrectly, based on the argument of the 
absence of I'i/'l^-term for D = 3 and of the absence of |-?/'pln|'?/^|-term for D = 2. 

More usual case of quenched disorder in Abelian-Higgs models has been considered in 
Refs. [75l [76] . This is the case when the parameter a in (1) has a random part, which 
obeys a Gaussian distribution. This disorder is of type random impurities, or, what is the 
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same, "random critical temperature [4]. In case of the so-called "uncorrelated", or, (6-) 
short range correlated quenched impurities, the system exhibits a spectacular competition 
between the impurities and gauge effects. The RG equations have a new stable FP of 
focal type [76] for dimensions D < 4, which exists for symmetry indices n > 1 and has 
a physical meaning at dimensions D > Dc{n) = 2{n + 36)/(n + 18) [76j. In the impure 
superconductor, the new focal FP occurs exactly in the domain of symmetry indices 
n < 365.9, where the HLM analysis [6] yields a lack of FP for the respective pure system. 
Having in mind the asymptotic nature of the e-expansions within the RG approach [U [77] , 
one may conclude that this focal FP governs the critical behavior at the real dimension 
D = 3 in conventional superconductors (n = 2), although the direct substitution of n = 2 
in the above expression for the "lower critical" dimension Dc{n) yields Dc(2) = 3.8. This 
problem has been further investigated in Ref. [78|. Long-range quenched disorder within 
the same RG approach to the model (1) has been considered in Ref. [79]. 

The effect of another relevant type of disorder, the so-called "random field disorder", on 
the phase transitions described by the GL functional (1) as been considered in Ref. [HO]. In 
both short-range and long-range random field correlations [H, [HH [82] a new stable FP has 
been obtained [80] . For short range random correlations this FP is stable below the upper 
critical dimension Du = 6 and for symmetry indices n > = 10. This FP describes 
a quite specific critical behavior. The situation for long range random correlations is 
more complicated and here we will advice the reader to follow the discussion in Ref. [80] . 
In general, the annealed disorder leaves the phase transition properties almost identical 
to those in the pure Abelian-Higgs model (1) whereas the quenched disorder gives a 
lower critical values of the symmetry index n, below which the weakly first order phase 
transition (HLM effect) exists. Thus one may conclude that the quenched disorder acts in 
a direction of "smearing" the fluctuation-driven first order phase transition, i.e. the HLM 
effect seems to be much weaker in such systems. 

3.4 Liquid crystals 

The first theoretical study of the HLM effect on the properties of the nematic-smectic 
A phase transition in liquid crystals has been performed in Refs. [22[ [23] . The main 
result of these investigations is the confirmation of the HLM effect in a close vicinity of 
this liquid crystal phase transition. Perhaps, because of the liquid crystal anisotropy, 
which explicitly enters in the propagator of the gauge field [23], here the critical value 
He = 38.17 of the symmetry index n, below which the weakly first order transition occurs, 
is lower than in pure superconductors (ric = 364.9) but it is still quite above the real value 
n = 2. Thus the RG predictions in Refs. [221 ESj in seventies of the previous century 
were in favor of a weakly first order phase transition whereas the experimentalists [S3l 
[M] at the same time claimed that the same nematic-smectic A phase transition is of 
second order (for triple and tricritical points in nematic-smectic- A-smectic-C systems, and 
the interrelationship to the present problem, see Ref. [85]), as well as later papers [86l 
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[STj. While in 3D superconductors the HLM effect is very weak and unobservable with 
available experimental techniques, the size of the same effect in 3D liquid crystals allows an 
observation in case of very precise experiment and elimination of obscuring effects as, for 
example, neighborhood to tricriticality and liquid crystal anisotropy. Experiments [HHl 
[89] consistent with the HLM effect has been carried out but, as noticed in Ref [90] . 
these experiments are at their resolution limits and the full implications of the HLM 
theory remains to be tested. On the other hand, experiments in Ref. [90] demonstrate a 
discrepancy with the simple approximation {ip =const) for the smectic-A filed within the 
mean-field like approach (Sec. II). 

3.5 Miscellaneous 

Quantum effects. The first account of quantum correlations ( "fluctuations" ) [H [9ll [92] 
on the properties of the phase transition to superconducting state has been performed in 
Ref. [93]. It has been shown [93] that the dynamical critical exponent z produced by the 
quantum correlations at finite temperatures (T > 0) is given by z = 2 + 18e/n. This is 
purely gauge result because, as shown in ref. [M], the critical dynamics of a superconductor 
in case of neglecting the local gauge effects {A = 0) is identical to that described by the 
time dependent GL model (TDGL) in case of lack of conservation laws [95] ; a prediction of 
the latter result has been given for the first time in Ref. [91] . Another study [96] of a gauge 
model of type (1) is intended to the study of the quantum phase transition (Tcq = 0) in 
granular superconductors. The simultaneous effect of the local U{1) symmetry, disorder 
and quantum fluctuations has been considered in the study [971 EH] of disordered electronic 
systems at T = 0. The problem was further extended to the treatment of quantum phase 
transitions in underdoped high-temperature superconductors [99 | llOOj . 

Complex systems. The RG investigation [101] of systems with superconducting and other 
(non- magnetic) orderings shows that the gauge field A{x) leads to a drastic modification of 
the critical behavior in close vicinity of bicritical and tetracritical points (for such points, 
see, e.g., Ref. [21 [1]). Near these multicritical points the superconducting fluctuations 
can be enhanced by the fluctuations of another ordering field and, hence, the HLM effect 
in such complex systems seems to be stronger than in the simpler case discussed so far 
[102] . The coupling between the superconducting field ip and the magnetization mode 
M{x) in models of ferromagnetic superconductors intended to describe the coexistence 
of superconductivity and ferromagnetism in ternary compounds has been investigated by 
RG in Ref. [1031 1104] . It has been shown again that the gauge field A{x) produces a 
weakly first order phase transition enhanced by the M-fluctuations. The same problem 
has been further opened in Ref. [105] by a lattice version of the GL model and dual argu- 
ments [ni] with the conclusion that under certain circumstances the phase transition from 
the disordered phase to the phase of coexistence of ferromagnetism and superconductivity 
should be, at least in type II superconductors, of second order in contrast to the RG pre- 
diction [1031 1104] . Another investigation of coupled magnetic and superconducting order 
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parameters, where the gauge effects are also relevant, is intended to describe a state of 
vortex solid in rear-earth and borocarbide compounds |106l I107j . 

Unconventional superconductors. Another interesting problem is the behavior of uncon- 
ventional superconductors [55], 1108] where the field ip is a. complex vector: ip = {ipa, a = 
2,3}, i.e., n = 4 or 6, depending on the type of the unconventional Cooper pairing. It 
has been shown |109j . that the HLM FP in this case is unstable towards the parameters, 
describing the anisotropy of the Cooper pair and the crystal anisotropy, and the new FP 
points that appear in the theory are unstable for usual values of the symmetry index 
(n = 4,6). It has been shown by analytical calculations [109] that one of the FP points 
that appear in the RG equations exists for n > Uc = 10988 and is stable for some values 
above Uc, at least, in the Hartree limit {n oo). Owing to these results, and mainly 
due to very high value of n^, a conclusion has been made that the HLM effect should be 
more pronounced in unconventional superconductors [109] . A next step along this path 
of RG studies has been made in Ref. |110j . where unconventional superconductors with 
quenched impurities have been considered. 

External field effect. In real experiments the external magnetic field can hardly be com- 
pletely eliminated. The study [38j of the HLM effect in the wider context of nonzero 
external magnetic field Hq shows that the weakly first phase transition discussed so far 
can be obtained in Landau gauge from a renormalized theory of the equation of state in 
one loop order; see also Ref. [lEj. The phase transition at the second critical magnetic 
field Hc2 has been studied by a e = (6 — rf) expansion within the one loop approximation 
with the conclusion that the magnetic fluctuation effects should produce a fluctuation 
induced first order phase transition [HI]. Further investigations of this problem have 
been performed in the Hartree limit {n 0) |112l I113j . and in higher orders of the loop 
expansion |114l 1115] : for comments, see Refs. |116|, 1117] . 
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